Characterization of the swede midge, Contarinia nasturtii, first instar larval salivary gland transcriptome.

Proteins in saliva of gall-forming insect larvae govern insect-host plant interactions. Contarinia nasturtii, the swede midge, is a pest of brassicaceous vegetables (cabbage, cauliflower, broccoli) and canola. We examined the salivary gland (SG) transcriptome of first instar larvae reared on Brassica napus and catalogued genes encoding secreted proteins that may contribute to the initial stages of larval establishment, the synthesis of plant growth hormones, extra-oral digestion and evasion of host defenses. A significant portion of the secreted proteins with unknown functions were unique to C. nasturtii and were often members of larger gene families organized in genomic clusters with conservation patterns suggesting that they are undergoing selection.

Horizontally transmitted parasitoid killing factor shapes insect defense to parasitoids.

Interkingdom competition occurs between hymenopteran parasitoids and insect viruses sharing the same insect hosts. It has been assumed that parasitoid larvae die with the death of the infected host or as result of competition for host resources. Here we describe a gene family, parasitoid killing factor (pkf), that encodes proteins toxic to parasitoids of the Microgastrinae group and determines parasitism success. Pkfs are found in several entomopathogenic DNA virus families and in some lepidopteran genomes. We provide evidence of equivalent and specific toxicity against endoparasites for PKFs found in entomopoxvirus, ascovirus, baculovirus, and Lepidoptera through a mechanism that elicits apoptosis in the cells of susceptible parasitoids. This highlights the evolutionary arms race between parasitoids, viruses, and their insect hosts.

De Novo Whole-Genome Assembly of the Swede Midge (Contarinia nasturtii), a Specialist of Brassicaceae, Using Linked-Read Sequencing.

The swede midge, Contarinia nasturtii, is a cecidomyiid fly that feeds specifically on plants within the Brassicaceae. Plants in this family employ a glucosinolate-myrosinase defense system, which can be highly toxic to nonspecialist feeders. Feeding by C. nasturtii larvae induces gall formation, which can cause substantial yield losses thus making it a significant agricultural pest. A lack of genomic resources, in particular a reference genome, has limited deciphering the mechanisms underlying glucosinolate tolerance in C. nasturtii, which is of particular importance for managing this species. Here, we present an annotated, scaffolded reference genome of C. nasturtii using linked-read sequencing from a single individual and explore systems involved in glucosinolate detoxification. The C. nasturtii genome is similar in size and annotation completeness to that of the Hessian fly, Mayetiola destructor, but has greater contiguity. Several genes encoding enzymes involved in glucosinolate detoxification in other insect pests, including myrosinases, sulfatases, and glutathione S-transferases, were found, suggesting that C. nasturtii has developed similar strategies for feeding on Brassicaceae. The C. nasturtii genome will, therefore, be integral to continued research on plant-insect interactions in this system and contribute to effective pest management strategies.

Examining population structure of a bertha armyworm, Mamestra configurata (Lepidoptera: Noctuidae), outbreak in western North America: Implications for gene flow and dispersal.

The bertha armyworm (BAW), Mamestra configurata, is a significant pest of canola (Brassica napus L. and B. rapa L.) in western North America that undergoes cyclical outbreaks every 6-8 years. During peak outbreaks millions of dollars are spent on insecticidal control and, even with control efforts, subsequent damage can result in losses worth millions of dollars. Despite the importance of this pest insect, information is lacking on the dispersal ability of BAW and the genetic variation of populations from across its geographic range which may underlie potential differences in their susceptibility to insecticides or pathogens. Here, we examined the genetic diversity of BAW populations during an outbreak across its geographic range in western North America. First, mitochondrial cytochrome oxidase 1 (CO1) barcode sequences were used to confirm species identification of insects captured in a network of pheromone traps across the range, followed by haplotype analyses. We then sequenced the BAW genome and used double-digest restriction site associated DNA sequencing, mapped to the genome, to identify 1000s of single nucleotide polymorphisms (SNP) markers. CO1 haplotype analysis identified 9 haplotypes distributed across 28 sample locations and three laboratory-reared colonies. Analysis of genotypic data from both the CO1 and SNP markers revealed little population structure across BAW's vast range. The CO1 haplotype pattern showed a star-like phylogeny which is often associated with species whose population abundance and range has recently expanded and combined with pheromone trap data, indicates the outbreak may have originated from a single focal point in central Saskatchewan. The relatively recent introduction of canola and rapid expansion of the canola growing region across western North America, combined with the cyclical outbreaks of BAW caused by precipitous population crashes, has likely selected for a genetically homogenous BAW population adapted to this crop.

Role of the peritrophic matrix in insect-pathogen interactions.

The peritrophic matrix (PM) is an acellular chitin and glycoprotein layer that lines the invertebrate midgut. The PM has long been considered a physical as well as a biochemical barrier, protecting the midgut epithelium from abrasive food particles, digestive enzymes and pathogens infectious per os. This short review will focus on the latter function, as a barrier to pathogens infectious per os. We focus on the evidence confirming the role of the PM as protective barrier against pathogenic microorganisms of insects, mainly bacteria and viruses, as well as the evolution of a variety of mechanisms used by pathogens to overcome the PM barrier.

Proteomics analysis of Trichoplusia ni midgut epithelial cell brush border membrane vesicles.

The insect midgut epithelium is composed of columnar, goblet, and regenerative cells. Columnar epithelial cells are the most abundant and have membrane protrusions that form the brush border membrane (BBM) on their apical side. These increase surface area available for the transport of nutrients, but also provide opportunities for interaction with xenobiotics such as pathogens, toxins and host plant allelochemicals. Recent improvements in proteomic and bioinformatics tools provided an opportunity to determine the proteome of the T. ni BBM in unprecedented detail. This study reports the identification of proteins from BBM vesicles (BBMVs) using single dimension polyacrylamide gel electrophoresis coupled with multi-dimensional protein identification technology. More than 3000 proteins were associated with the BBMV, of which 697 were predicted to possess either a signal peptide, at least one transmembrane domain or a GPI-anchor signal. Of these, bioinformatics analysis and manual curation predicted that 185 may be associated with the BBMV or epithelial cell plasma membrane. These are discussed with respect to their predicted functions, namely digestion, nutrient uptake, cell signaling, development, cell-cell interactions, and other functions. We believe this to be the most detailed proteomic analysis of the lepidopteran midgut epithelium membrane to date, which will provide information to better understand the biochemical, physiological and pathological processes taking place in the larval midgut.

The complete genome sequence of a second alphabaculovirus from the true armyworm, Mythimna unipuncta: implications for baculovirus phylogeny and host specificity.

The Mythimna unipuncta nucleopolyhedrovirus isolate KY310 (MyunNPV-KY310) is an alphabaculovirus isolated from a true armyworm (Mythimna unipuncta) population in Kentucky, USA. Occlusion bodies of this virus were examined by electron microscopy and the genome sequence was determined by 454 pyrosequencing. MyunNPV-KY310 occlusion bodies consisted of irregular polyhedra measuring 0.8-1.8 µm in diameter and containing multiple virions, with one to six nucleocapsids per virion. The genome sequence was determined to be 156,647 bp with a nucleotide distribution of 43.9% G+C. 152 ORFs and six homologous repeat (hr) regions were annotated for the sequence, including the 38 core genes of family Baculoviridae and an additional group of 26 conserved alphabaculovirus genes. BLAST queries and phylogenetic inference confirmed that MyunNPV-KY310 is most closely related to the alphabaculovirus Leucania separata nucleopolyhedrovirus isolate AH1, which infects Mythimna separata. In contrast, MyunNPV-KY310 did not exhibit a close relationship with Mythimna unipuncta nucleopolyhedrovirus isolate #7, an alphabaculovirus from the same host species. MyunNPV-KY310 lacks the gp64 envelope glycoprotein, which is a characteristic of group II alphabaculoviruses. However, this virus and five other alphabaculoviruses lacking gp64 are placed outside the group I and group II clades in core gene phylogenies, further demonstrating that viruses of genus Alphabaculovirus do not occur in two monophyletic clades. Potential instances of MyunNPV-KY310 ORFs arising by horizontal transfer were detected. Although there are now genome sequences of four different baculoviruses from M. unipuncta, comparison of their genome sequences provides little insight into the genetic basis for their host specificity.

The Effect of Diet on Midgut and Resulting Changes in Infectiousness of AcMNPV Baculovirus in the Cabbage Looper, Trichoplusia ni.

Insecticide resistance has been reported in many important agricultural pests, and alternative management methods are required. Baculoviruses qualify as an effective, yet environmentally benign, biocontrol agent but their efficacy against generalist herbivores may be influenced by diet. However, few studies have investigated the tritrophic interactions of plant, pest, and pathogen from both a gene expression and physiological perspective. Here we use microscopy and transcriptomics to examine how diet affects the structure of peritrophic matrix (PM) in Trichoplusia ni larvae and consequently their susceptibility to the baculovirus, AcMNPV. Larvae raised on potato leaves had lower transcript levels for chitinase and chitin deacetylase genes, and possessed a thicker and more multi-layered PM than those raised on cabbage or artificial diet, which could contribute to their significantly lower susceptibility to the baculovirus. The consequences of these changes underline the importance of considering dietary influences on pathogen susceptibility in pest management strategies.

The complete genome sequence of a third distinct baculovirus isolated from the true armyworm, Mythimna unipuncta, contains two copies of the lef-7 gene.

A baculovirus isolate from a USDA Forest Service collection was characterized by electron microscopy and analysis of its genome sequence. The isolate, formerly referred to as Pseudoletia (Mythimna) sp. nucleopolyhedrovirus #7 (MyspNPV#7), was determined by barcoding PCR to derive from the host species Mythimna unipuncta (true armyworm) and was renamed Mythimna unipuncta nucleopolyhedrovirus #7 (MyunNPV#7). The occlusion bodies (OBs) and virions exhibited a size and morphology typical for OBs produced by the species of genus Alphabaculovirus, with occlusion-derived virions consisting of 2-5 nucleocapsids within a single envelope. The MyunNPV#7 genome was determined to be 148,482 bp with a 48.58% G+C nucleotide distribution. A total of 159 ORFs of 150 bp or larger were annotated in the genome sequence, including the 38 core genes of family Baculoviridae. The genome contained six homologous repeat regions (hrs) consisting of multiple copies of a 34-bp imperfect palindrome. Phylogenetic inference from concatenated baculovirus core gene amino acid sequence alignments placed MyunNPV#7 with group II alphabaculoviruses isolated from other armyworm and cutworm host species of lepidopteran family Noctuidae. MyunNPV#7 could be distinguished from other viruses in this group on the basis of differences in gene content and order. Pairwise nucleotide distances suggested that MyunNPV#7 represents a distinct species in Alphabaculovirus. The MyunNPV#7 genome was found to contain two copies of the late expression factor-7 (lef-7) gene, a feature not reported for any other baculovirus genome to date. Both copies of lef-7 encoded an F-box domain, which is required for the function of LEF-7 in baculovirus DNA replication.

Baculoviruses require an intact ODV entry-complex to resist proteolytic degradation of per os infectivity factors by co-occluded proteases from the larval host.

Baculoviruses orally infect caterpillars in the form of occlusion-derived viruses (ODVs). The ODV-envelope contains a number of proteins which are essential for oral infectivity, called per os infectivity factors (PIFs). Most of these PIFs are involved in the formation of an ODV-entry complex that consists of a stable core, formed by PIF1, PIF2, PIF3 and PIF4, and the more loosely associated PIFs P74 (PIF0) and P95 (PIF8). PIF1, PIF2 and PIF3 are essential for formation of the stable core, whereas deletion of the pif4 gene results in the formation of a smaller complex. P74 is not needed for formation of the stable core. We show here in larva-derived ODVs of the Autographa californica multicapsid nucleopolyhedrovirus that PIF-proteins are degraded by host-derived proteases after deletion of a single pif-gene. Constituents of the stable core-complex appeared to be more resistant to proteases as part of the complex than as monomer, as in ODVs of a p74 deletion mutant only the stable core was found but no PIF monomers. When the stable core lacks PIF4, it lost its proteolytic resistance as the resulting smaller core complex was degraded in a pif4 deletion mutant. We also identified PIF6 as a loosely associated component of the entry complex that appeared nevertheless important for the proteolytic resistance of the stable core, which was degraded after deletion of pif6. We conclude from these results that an intact entry-complex in the ODV-envelope is prerequisite for proteolytic resistance of PIF-proteins under the alkaline conditions of the larval midgut.

Lacanobia oleracea nucleopolyhedrovirus (LaolNPV): A new European species of alphabaculovirus with a narrow host range.

During an insect sampling program in alfalfa crops near Montpellier, France in 2011, Lacanobia oleracea larvae were collected that died due to nucleopolyhedrovirus infection (LaolNPV). This virus was subjected to molecular and biological characterization. The virus was a multiple nucleocapsid NPV that showed similar restriction profiles to Mamestra configurata NPV-A (MacoNPV-A) but with significant differences. Polypeptide analysis demonstrated similar proteins in occlusion bodies and occlusion derived virions, to those observed in NPVs from Mamestra spp. Terminal sequencing revealed that the genome organization shared similarity with that of MacoNPV-A. The most homologous virus was MacoNPV-A 90/2 isolate (95.63% identity and 96.47% similarity), followed by MacoNPV-A 90/4 strain (95.37% and 96.26%), MacoNPV-B (89.21% and 93.53%) and M. brassicae MNPV (89.42% and 93.74%). Phylogenetic analysis performed with lef-8, lef-9, polh and a concatenated set of genes showed that LaolNPV and the Mamestra spp. NPVs clustered together with HaMNPV, but with a closer genetic distance to MacoNPV-A strains. The Kimura 2-parameter (K-2-P) distances of the complete genes were greater than 0.05 between LaolNPV and the MbMNPV/MacoNPV-B/HaMNPV complex, which indicates that LaolNPV is a distinct species. K-2-P distances were in the range 0.015-0.050 for comparisons of LaolNPV with MacoNPV-A strains, such that additional biological characteristics should be evaluated to determine species status. While MacoNPV-A was pathogenic to seven lepidopteran species tested, LaolNPV was only pathogenic to Chrysodeixis chalcites. Given these findings, Lacanobia oleracea nucleopolyhedrovirus should be considered as a new species in the Alphabaculovirus genus.

The Complete Genome Sequence of a Second Distinct Betabaculovirus from the True Armyworm, Mythimna unipuncta.

The betabaculovirus originally called Pseudaletia (Mythimna) sp. granulovirus #8 (MyspGV#8) was examined by electron microscopy, host barcoding PCR, and determination of the nucleotide sequence of its genome. Scanning and transmission electron microscopy revealed that the occlusion bodies of MyspGV#8 possessed the characteristic size range and morphology of betabaculovirus granules. Barcoding PCR using cytochrome oxidase I primers with DNA from the MyspGV#8 collection sample confirmed that it had been isolated from the true armyworm, Mythimna unipuncta (Lepidoptera: Noctuidae) and therefore was renamed MyunGV#8. The MyunGV#8 genome was found to be 144,673 bp in size with a nucleotide distribution of 49.9% G+C, which was significantly smaller and more GC-rich than the genome of Pseudaletia unipuncta granulovirus H (PsunGV-H), another M. unipuncta betabaculovirus. A phylogeny based on concatenated baculovirus core gene amino acid sequence alignments placed MyunGV#8 in clade a of genus Betabaculovirus. Kimura-2-parameter nucleotide distances suggested that MyunGV#8 represents a virus species different and distinct from other species of Betabaculovirus. Among the 153 ORFs annotated in the MyunGV#8 genome, four ORFs appeared to have been obtained from or donated to the alphabaculovirus lineage represented by Leucania separata nucleopolyhedrovirus AH1 (LeseNPV-AH1) during co-infection of Mythimna sp. larvae. A set of 33 ORFs was identified that appears only in other clade a betabaculovirus isolates. This clade a-specific set includes an ORF that encodes a polypeptide sequence containing a CIDE_N domain, which is found in caspase-activated DNAse/DNA fragmentation factor (CAD/DFF) proteins. CAD/DFF proteins are involved in digesting DNA during apoptosis.

Autographa californica Multiple Nucleopolyhedrovirus AC83 is a Per Os Infectivity Factor (PIF) Protein Required for Occlusion-Derived Virus (ODV) and Budded Virus Nucleocapsid Assembly as well as Assembly of the PIF Complex in ODV Envelopes.

Baculovirus occlusion-derived virus (ODV) initiates infection of lepidopteran larval hosts by binding to the midgut epithelia, which is mediated by per os infectivity factors (PIFs). Autographa californica multiple nucleopolyhedrovirus (AcMNPV) encodes seven PIF proteins, of which PIF1 to PIF4 form a core complex in ODV envelopes to which PIF0 and PIF6 loosely associate. Deletion of any pif gene results in ODV being unable to bind or enter midgut cells. AC83 also associates with the PIF complex, and this study further analyzed its role in oral infectivity to determine if it is a PIF protein. It had been proposed that AC83 possesses a chitin binding domain that enables transit through the peritrophic matrix; however, no chitin binding activity has ever been demonstrated. AC83 has been reported to be found only in the ODV envelopes, but in contrast, the Orgyia pseudotsugata MNPV AC83 homolog is associated with both ODV nucleocapsids and envelopes. In addition, unlike known pif genes, deletion of ac83 eliminates nucleocapsid formation. We propose a new model for AC83 function and show AC83 is associated with both ODV nucleocapsids and envelopes. We also further define the domain required for nucleocapsid assembly. The cysteine-rich region of AC83 is also shown not to be a chitin binding domain but a zinc finger domain required for the recruitment or assembly of the PIF complex to ODV envelopes. As such, AC83 has all the properties of a PIF protein and should be considered PIF8. In addition, pif7 (ac110) is reported as the 38th baculovirus core gene.IMPORTANCE ODV is essential for the per os infectivity of the baculovirus AcMNPV. To initiate infection, ODV binds to microvilli of lepidopteran midgut cells, a process which requires a group of seven virion envelope proteins called PIFs. In this study, we reexamined the function of AC83, a protein that copurifies with the ODV PIFs, to determine its role in the oral infection process. A zinc finger domain was identified and a new model for AC83 function was proposed. In contrast to previous studies, AC83 was found to be physically located in both the envelope and nucleocapsid of ODV. By deletion analysis, the AC83 domain required for nucleocapsid assembly was more finely delineated. We show that AC83 is required for PIF complex formation and conclude that it is a true per os infectivity factor and should be called PIF8.

Microscopic investigation of AcMNPV infection in the Trichoplusia ni midgut.

Autographa californica multiple nucleopolyhedrovirus (AcMNPV) is the type species for the genus Alphabaculovirus in the family Baculoviridae. In nature, AcMNPV infection begins with ingestion of viral occlusion bodies (OBs) from which occlusion-derived viruses (ODV) are released to infect midgut cells. This study explored the early stages of Trichoplusia ni midgut infection using recombinant viruses expressing green fluorescent protein (GFP) and/or a VP39-mCherry fusion protein under the control of early and late promoters, respectively. Using a recombinant ie1:GFP virus, the anterior midgut region was identified as the predominant site for primary infection. Infection of midguts using the GFP-VP39mCherry-dual labelled recombinant virus revealed that active viral replication and cell-to-cell spread was required for the formation of infection foci and the subsequent spread to uninfected midgut cells and tracheoblasts. The spread of the infection from primary infected cells to secondary cells within the midgut was shown to be dependent upon the membrane fusion protein GP64.

MacoNPV baculovirus midgut-specific gene expression during infection of the bertha armyworm, Mamestra configurata.

Baculoviruses have two forms, occlusion derived virus (ODV) which is responsible for primary infection in host midgut tissue and budded virus (BV), which infects all other host tissues during secondary infection. This study examined the primary infection by ODV of midgut cells of bertha armyworm Mamestra configurata fourth instar larvae and measured the expression of viral genes over a time course of infection. Both digital PCR and RNA sequencing methods showed the profile of transcription to be different from those produced by AcMNPV BV infection of in vitro cell cultures. This included having unique collections of genes expressed early, as well as much greater late gene expression of p6.9 and much reduced expression of polh and p10. These differences likely reflect characteristics unique to the critical step of in vivo midgut cell infection, and provide insights into the processes that regulate viral gene expression in different host tissues.

Genome Sequence of an Alphabaculovirus Isolated from the Oak Looper, Lambdina fiscellaria, Contains a Putative 2-Kilobase-Pair Transposable Element Encoding a Transposase and a FLYWCH Domain-Containing Protein.

The genome sequence of an alphabaculovirus isolated from Lambdina fiscellaria indicated that it is a novel member of a group II lineage. A putative transposable element was identified that contained two genes, including a transposase ortholog. These genes were most closely related to genes of the pea aphid, Acyrthosiphon pisum.

A Distinct Group II Alphabaculovirus Isolated from a Peridroma Species.

The genome sequence of an alphabaculovirus isolated from a Peridroma species indicated that it is a novel member of a group II lineage most closely related to alphabaculoviruses from Spodoptera exigua and Agrotis segetum. It contains a genome of 151,110 nucleotides (nt), with a G+C content of 53.3%.

Genome Sequence of an Alphabaculovirus Isolated from Choristoneura murinana.

The genome sequence of a baculovirus from Choristoneura murinana is 124,689 bp, with a G+C content of 50%, and contains 148 putative open reading frames. The virus is a member of the group I alphabaculoviruses and is most closely related to several other viruses that infect Choristoneura species.

Mamestra configurata nucleopolyhedrovirus-A transcriptome from infected host midgut.

Infection of an insect by a baculovirus occurs in two distinct phases, an initial infection of host midgut by occlusion-derived virions (ODVs) and subsequent systemic infection of other tissues by budded virions (BV). A vast majority of investigations of the infection process have been restricted to cell culture studies using BV that emulate the systemic phase of infection. This is one of the first studies to investigate baculovirus gene expression in ODV infected midgut cells. We have focused on the critical first phase of in vivo infection by Mamestra configurata nucleopolyhedrovirus-A in M. configurata larvae, using qPCR and RNAseq mass sequencing to measure virus gene expression in midgut cells. The earliest genes detected by each method had significant overlap, including known early genes as well as genes unique to MacoNPV-A and genes of unknown function. The RNAseq data also revealed a large range of expression levels across all ORFs, which could not be measured using qPCR. This dataset provides a first whole genome transcriptomic analysis of viral genes required for virus infection in vivo and will provide the basis for functionally analyzing specific genes that may be critical elements in baculovirus midgut infectivity and host range.

The genome of a baculovirus isolated from Hemileuca sp. encodes a serpin ortholog.

The genome sequence of a baculovirus from Hemileuca sp. was determined. The genome is 140,633 kb, has a G+C content of 38.1 %, and encodes 137 putative open-reading frames over 50 amino acids. 126 of these ORFs showed similarity to other baculovirus genes in the database including all 37 core genes. Of the remaining 11 predicted genes, one is related to a lepidopteran serpin gene. This is the first report of a baculovirus encoding a member of this family of serine protease inhibitors, and to our knowledge the first report of a viral serpin outside the Poxviridae. The genome also contained three homologous repeat sequences. Phylogenetic analysis indicated that the virus is a group II Alphabaculovirus and belongs to a lineage that includes Orgyia leucostigma, Ectropis obliqua, Apocheima cinerarium, and Euproctis pseudoconspersa nucleopolyhedroviruses.